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An experimental study of Taylor's hypothesis in a quasi-two-dimensional turbulent soap 
film is presented. A two probe laser Doppler velocimeter enables a non-intrusive simultane- 
ous measurement of the velocity at spatially separated points. The breakdown of Taylor's 
hypothesis is quantified using the cross correlation between two points displaced in both 
space and time; correlation is better than 90% for scales less than the integral scale. A 
quantitative study of the decorrelation beyond the integral scale is presented, including an 
analysis of the failure of Taylor's hypothesis using techniques from predictability studies 
of turbulent flows. Our results are compared with similar studies of 3D turbulence. 
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I. INTRODUCTION 

In a 1938 paper on the statistics of turbulence, G. I. Taylor presented an assumption 
from which he could infer the spatial structure of a turbulent velocity field from a single 
point measurement of its temporal fluctuation |1|. This assumption, known as Taylor's 
hypothesis or the frozen turbulence assumption, relies on the existence of a large mean flow 
which translates the fluctuations past the stationary probe in a time short compared to 
the evolution time of the turbulence. The experimental measurements treated by Taylor 
were made on the turbulence generated behind a stationary grid in a wind tunnel, and his 
hypothesis has become a standard technique employed in similar experiments which inform 
our current views on turbulence (see for example Refs. 0, [Q, and [Q). The importance of 
this hypothesis stems from the fact that most turbulence theories are framed in terms of the 
spatial structure of the velocity field @J^. 

In practical terms, the limits of Taylor's hypothesis are determined by how large the 
mean velocity must be relative to the fluctuations. Recently Yakhot has pointed out 
that the corrections to Taylor's hypothesis could well be of the same order as corrections 
to the standard model of turbulence, Kolmogorov's 1941 theory There is at present, 

however, no firm theoretical derivation of the hypothesis, and thus no fashion to evaluate 
its reliability, calculate higher order corrections, or predict in what way it will break down. 
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A few theoretical discussions do exist P-p!T[], and a promising direction has recently been 



taken by Hayot & Jayaprakash for the Burgers equation [|12|. The treatment by Lumley 



in particular offers corrections to statistical measures of spatial gradients for non-negligible 



velocity fluctuations [T^ . Still a general theoretical framework is lacking. 

Experimentally there have been many studies of Taylor's hypothesis, all of which have 
treated three dimensional (3D) turbulence; we give a brief overview in Section II. Our 
experiments, however, are performed on the quasi-two dimensional flow of a soap film. Two- 
dimensional (2D) fluid flows occur in many physical situations, mostly due to the effects of 
rotation or stratification in the atmosphere and ocean [T^|. Turbulence in 2D is different 



from 3D in several ways, largely due to the absence of vortex stretching in 2D [|14|. This 
means that the squared vorticity (called enstrophy) becomes a nearly-conserved quantity 
in 2D, like the energy, and thus two cascades are expected: a direct cascade of enstrophy 
to smaller scales, and an inverse cascade of energy to larger scales [0; for decaying 2D 
turbulence the inverse cascade is apparently absent |n,rB[. Although Taylor's hypothesis is 
also an important assumption in the study of 2D turbulence, to our knowledge it has never 
been tested, nor is it clear that the hypothesis should do relatively better or worse than in 
3D. 

It may be helpful at this point to discuss the essential differences between 2D and 3D 
turbulence in order to better relate the present measurements to prior three dimensional 
tests of the Taylor hypothesis. In three dimensions the vorticity vector uj{x, y, z, t) can point 
in any direction, but in 2D it is restricted to be perpendicular to the x, y plane of the flow. 
This fact alone assures that in inviscid flows, the enstrophy f2 = | < cj^ > is a constant of 
the motion, in addition to kinetic energy conservation K = ^ < > (the angular brackets 
designate an appropriate average). From the Navier-Stokes equation it follows that in 2D, 
vorticity cannot be amplified (or attenuated except by viscous damping) by a velocity field 
gradient. The existence of vortex stretching in 3D is intimately related to the energy cascade 
from large scales to small, with the process controlled by the rate e at which kinetic energy 
is injected at large scales. 

In 2D, one expects that energy injected at intermediate scale, rinj will be transferred to 
larger scales [0,0 and dissipated at the boundaries of the system. This inverse cascade 
process is not expected to be local in k-space. For scales r < rj„j, the small-scale velocity 
fluctuations < 15^(^)1 >, which carry little of the energy, are expected to cascade down to 
the dissipative scale under the control of the enstrophy injection rate f3 = dQ/dt (the energy 
injection rate e plays the corresponding role in the 2D inverse cascade). 

It seems reasonable that the absence of vortex stretching in 2D will increase the likelihood 
of a velocity fluctuation being transported intact to a distant downstream point, thereby 
enhancing the validity of Taylor's hypothesis over 3D turbulence. Taylor's hypothesis can 
also fail when a local fluctuation is transported laterally into the flow path between adjacent 
points separated by the distance Ax = UoAt. For purely geometrical reasons, this would 
seem to be a less likely occurence in 2D than in 3D, so that the frozen turbulence relation 
6v{Ax) = 6v{UQt) should presumably hold out to larger values of Ax. It is harder for us 
to assess the impact of 2D nonlocality on the validity of the Taylor hypothesis in 2D as 
compared to 3D. 

Our experiment uses a laser Doppler velocimeter (LDV) with two probes to nonintrusively 
examine Taylor's hypothesis in a turbulent soap film. We differentiate between two different 
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aspects of the hypothesis. By Taylor's hypothesis of coherence we mean the assertion that 
the velocity field is unchanged as it is advected downstream: that v{x,t) is identical to 
v{x + Ax, t + To), where tq = Ax/Uq, with Uq being the mean velocity and Ax the distance 
between the two points . For a clear illustration of this, glance ahead to Figure This 
coherence hypothesis is unquestionably an approximation and must fail as Aa; becomes large; 
here we quantify this failure, and relate it to the predictability problem for turbulent flows. 
We measure the breakdown of the coherence hypothesis via the correlation between two 
points in the flow, displaced in both position and time [p!^-p3|. We flnd nevertheless that 
the expectation value of the lowest six moments of the longitudinal velocity difference is 
unaffected by the loss of coherence. Thus time correlation statistics of velocity fluctuations 
appear to be the same as spatial ones; we will refer to this as Taylor's statistical hypothesis, 
which is implied by the coherence hypothesis, but does not actually require it. 



II. PREVIOUS STUDIES OF TAYLOR'S HYPOTHESIS 

Taylor's hypothesis has been exposed to many experimental tests in three dimensional 
turbulence, and we do not intend to give an exhaustive survey here. These experimental 
studies can be divided into two broad categories, concerned with either correlations over 



flnite distances [p!^-p3|, or local spatial derivatives used in turbulent dissipation estimates 
P^^7|. It has long been appreciated that the validity of Taylor's frozen turbulence as- 
sumption requires the smallness of the turbulent intensity It, deflned as the ratio of rms 
velocity fluctuations to the mean flow speed Uq. Additionally, the mean shear rate and 
the viscous damping must be small in the range of spatial scales Ax being probed. In this 
section we give a sampling of the sort of work which has been done in 3D; a more detailed 
discussion of studies with which we directly compare our results is given in Section V. The 



reader is referred to introductory reviews in two recent articles p6| , p7| , which are somewhat 
complementary to what is given here. 

To test the application of Taylor's hypothesis to velocity correlations over flnite dis- 
tances, one approach is to compare measurements made at a single observation point with 
measurements made at points displaced downstream. The flrst such tests were made in 
a wind tunnel by Favre, Gaviglio, & Dumas |jl8|,|l9|. The measurements were performed 



within the turbulent boundary layer of a plate at various It up to 15%. They found that 
Taylor's statistical hypothesis is valid for measurements of the velocity correlation function 
R{Ax,t) = {vi{xi,t)v2{xi + Ax,t — r)) made not too close to the plate. Fisher & Davies 
pO| made careful measurements of the velocity correlation function R{Ax, r) in a jet. They 
found that the relation r = Ax/Uq was not well satisfled in the mixing region, where It is 
typically ~ 20%. These authors observed, as have many others, that the functional form of 
R changes with increasing Ax, and that this function is not very sharply peaked, as it would 
be if Taylor's coherence hypothesis were satisfled. Comte-Bellot & Corrsin measured 



R{Ax, t) for grid-generated turbulence in a wind tunnel, where both It and the mean shear 
rates are rather small. Though downstream decay causes the maximum value of R vs. r to 
decrease with increasing Ax, the correlation functions could still be collapsed onto the same 
functional form. 

One of the fundamental effects of turbulence is its enhancement of dissipation, and mea- 
surements of this require knowledge of spatial derivatives. Using Taylor's hypothesis allows 
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the time deriative of a single point measurement to be related to spatial derivatives, the 
simplest relation being dip/dx = Uodcp/dt, where could be a passive scalar concentration, 
temperature, a velocity component, or a product of velocity components. Kailasneth, Sreeni- 
vasan, & Saylor p5| studied a variety of turbulent systems and tested Taylor's hypothesis 
using a fluorescent dye in a jet, the heated wake of a cylinder, and the atmospheric boundary 
layer. They were interested in the probability density function of these scalars, and found 
that Taylor's statistical hypothesis worked well for conditional probability densities of the 
scalar fluctuations. Mi & Antonia ||2^ used a heated jet of air, with a turbulent intensity of 
about 26 %, to verify different theoretical relations between spatial and temporal derivatives 
of temperature, corrected for finite turbulent intensities. Dahm & Southerland compared 
2D spatial and spatio-temporal gradient fields of fluorescent dye in a water jet using a fast 
photodiode array. They found that Taylor's coherence hypothesis was only approximately 
verified for these fields. Piomelli, Balint, & Wallace p4| studied Taylor's hypothesis for 
various velocity derivatives and compared hot wire measurements, large eddy simulations, 
and direct numerical simulations of the Navier-Stokes equation for wall-bounded flows. Tay- 
lor's hypothesis was found to be in accord with the calculations and measurements made 
sufficiently far from the wall, where the mean shear is not excessive. In our experiments we 
have not treated the application of Taylor's hypothesis to gradients. 

All hot wire measurements are sufficiently intrusive that one must compensate or oth- 
erwise adjust for the perturbations produced by the wake of the upstream probe on the 
velocity measured at the downstream probe. Cenedese et al. p3[ avoided this problem by 
making velocity measurements with a Laser Doppler system. Only for small Ax did their 
correlation measurements satisfy Taylor's hypothesis very well, though admittedly It was 
rather high (13%). We discuss their results in more detail in Section V. 

In the present experiment. Laser Doppler velocimetry is also used to measure i?(Ax, r), so 
that there is no perturbation on the downstream probe. In contrast to the studies discussed 
above, however, we examine Taylor's hypothesis in a quasi-2D system, a flowing soap film; 
we have also studied the effects on the higher order moments of velocity differences. Before 
presenting our results, we discuss in some detail our experimental system. 



III. EXPERIMENTAL SETUP 

The use of soap films as convenient systems for the experimental study of 2D hydrody- 
namics began with the pioneering work of Couder and coworkers pB|-pn| and Gharib and 
Derango Our measurements are performed using a fiowing soap film apparatus devel- 
oped at the University of Pittsburgh by Kellay, Wu, & Goldburg [P2[, and Rutgers, Wu, & 



Goldburg ||33| , |34| ; we are using the latest version of this system, built by Rutgers. In our 
setup, a thin soap film several /im thick is allowed to fall between two taut plastic wires 
from an upper reservoir into a lower one, see Figure 0. The channel width is W = 6.2 cm 
over a distance of 120 cm, where the measurements are performed. Quasi-two-dimensional 
turbulence is generated behind a comb (tooth diameter 1 mm and spacing M = 3.8 mm) 
which perforates the film at a fixed height. The typical transit time between comb and lower 
reservoir is the order of a second. 

The turbulence generated in such a soap film decays downstream from the grid, exhibiting 
many aspects which agree with theories of 2D turbulence [^^-^,^,3B|, though there are also 
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some differences. It tlierefore seems wortliwliile to evaluate tlie experimental situation to 
date. Because the lateral dimensions of soap films are many orders of magnitude larger than 
their thickness, there would seem to be no doubt that the vorticity can indeed be regarded 
as a scalar quantity, so that vortex stretching is absent This would be a key requirement 
for the occurence two-dimensional turbulence in soap films. In addition, to be able to 
compare with theoretical and numerical results concerning incompressible turbulence, the 
two-dimensional compressibility of the film should be zero. This condition is clearly not fully 
met, since fluctuations in film thickness are visible to the eye through optical interference 
of light reflected from the front and back faces of the film |3D,JI,3B|. There are, however, 
recent experiments in which the two dimensional divergence D2 = V2v(a;, y) was measured 
by particle imaging velocimetery |Q. In a set-up very similar to that used here, D2 was 
measured to be 10-15% of the rms vorticity near the comb. Note also that because the 
velocity of peristaltic waves in a soap film is orders of magnitude larger than the turbulent 



velocity fluctuations BO], it is expected that the film may be regarded as incompressible 



from the point of view of this study. 

Another mitigating factor to the two-dimensionality of soap film flow is the friction 
between the film and the surrounding air. In a set of experiments in which the film was 
placed in a partial vacuum where the air pressure 3 % of the atmospheric value, it was found 
that the energy spectrum E{k) decayed for a decade in k as E{k) oc k~'^, with ( = 3.3± 0.3. 
This exponent had the same value in both the partial vacuum and at atmospheric pressure. 
The main effect of the reduced pressure was to magnify the magnitude of E{k) near the 



comb |35[. These measurements also showed that the total kinetic energy initially decayed 



downstream, but ultimately leveled off at a sufficient distance below the comb, which is 



expected theoretically at very high Reynolds numbers |]T6|. If the levelling off distance is 
called X* and the corresponding time t* = x*/Uq, then all measurements reported in this 
paper were made at values of t < t*, i.e. at distances from the comb too close for the levelling 
off to have occurred. 

In most (but not all) measurements of decaying turbulence in a soap film, only the en- 
strophy cascade is observed, i.e. the inverse energy cascade {E{k) ~ k~^^^) is not. Recently 
an experiment was performed in which turbulence was forced by an array of teeth parallel 
to the direction of flow There one finds evidence of both the inverse cascades 
and the inverse cascade spectrum, where k~^ . 



We use a commercial LDV system ||33] to measure the film velocity fluctuations ||38| , |39 



The soap solution (water and 2% commercial detergent by volume) is seeded with 1 /im 
polystyrene spheres at a volume fraction of about 10~^, and the data rate ranges from 1 
to 8 kHz. At a distance x = 8 cm below the comb, the mean and RMS velocities were 
typically Uq = 180 cm/s and vrms = (v'"^)^/'^ = 24 cm/s in the longitudinal (streamwise) 
direction, where v' = v — Uq. The turbulent intensity in these experiments was It — 0.14. 
This is the quantity which is assumed in Taylor's hypothesis to be small [|1],3, and we have 
explicitly chosen for our study a value of It which is not very small. The Reynolds number 
for the channel is Rew = UqW/u ~ 11,000, and for the comb Rcm = UqM/u ~ 700. The 
viscosity of a flowing soap film is not a well established quantity; here we use z/ = 0.1 cm^/s 
as measured in a 2D Couette viscometer by Martin & Wu Deviations from two- 



dimensionality caused by air friction appear not to affect the turbulence for the scales 
of interest here ^5 . 
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In order to test Taylor's hypothesis, two LDV heads are used at spatially separated 
points. Figure 1 shows the arrangement, with the downward (flow) direction defined as the 
X direction. One head (labeled "LDV 1") is kept fixed at Xi = 8 cm below the comb, while 
a second head ("LDV 2") is placed at various points ranging from 0:2 = 4 to 30 cm below 
the comb. The two probes are arranged to be directly in line with each other, so that the 
lower probe measures the same part of the fiow as the upper one, with a delay given by the 
transit time between them. 

Because the LDV only measures velocity when there is a scatterer in its measuring 
volume, the two probes do not in general measure velocity simultaneously. Therefore some 
binning or "simultaneity window" At is needed to perform statistical comparisons: two 
measurements are treated as simultaneous if they occur within At of one another. Here we 
use binning windows from At = 25 //s to 200 /xs, corresponding to frequencies 1/At from 
5 to 40 kHz, which are higher than the largest observed frequency in the velocity power 
spectrum. All measurements reported here are insensitive to small changes in At. 

One of the difficulties in examining the validity of Taylor's hypothesis in decaying tur- 
bulence is that the "small parameter" It is not constant, but decreases downstream as the 
turbulence decays. In Fig. |]a we plot our turbulent intensity as a function of distance from 
the comb. Thus, although we measure the correlation of velocity fiuctuations relative to 
Xi = 8 cm, where It — 0.14, the actual turbulent intensity affecting the velocity field down- 
stream is always less than 0.14. This would seem to significantly complicate the matter. 

How does the turbulence decay in our system? A standard result from 3D decaying grid 
turbulence is that the inverse square of the turbulent intensity depends on the distance from 
the grid as: /^"^ = A(x/M — B)^, where the dimensionless constants are typically found to 



be A ~ 130 - 150 and 5 ~ 3 - 20 for /? = 1 |41|, or A ~ 20 and B ~ 3.5 for /? = 1.25 [22 



note that B is the effective position of the origin for this scaling in units of M. Assuming 
that P = 1.25 means that If^'^ should be a linear function of x/M; however we find that by 
taking I^^'^ we get the best linear plot vs. x/M, as shown in Fig. The line corresponds 
to 

1 / X 



with i? = for our fit, which means that the virtual origin is located at the position of the 
grid itself. Note that we do not measure far enough downstream |^ to see any evidence of 
the kinetic energy saturation as the turbulence decays downstream {It ~ constant). 



IV. RESULTS 

A. Testing Taylor's Coherence Hypothesis 

Figure P shows the velocity fluctuations measured by two probes with Ax = 0.5 cm. For 
this small separation, Taylor's hypothesis is clearly a good estimate: the two velocity traces 
are nearly identical except for a small shift in time, which should correspond to the transit 
time across the spatial separation Ax. To test whether the velocity trace is translated spa- 
tially without evolving dynamically, we measure the cross-correlation Ci2(t, Ax, Xi) between 
the two probes: 
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, _ {vi{xi,t)v2{xi + Ax,t- t)) 
Ci2(r, Ax,xi) = . (1) 

VlRMS X V2RMS 

Here f t) and V2{x2, t) are the two velocities measured by the probes, virms are the RMS 
velocity fluctuations, Ax = X2 — Xi, and the brackets () denote a time average. Note that 
Ci2 is a function not only of delay time r and separation Ax, but of the absolute location of 
the flrst probe Xi. This last dependence comes from the fact that the turbulence is decaying. 
In this study we flx xi = 8 cm, and ignore this dependence. 



Figure 10 shows Ci2(t, Ax) vs. r for several different separations Ax. As expected, there 



is a well-deflned maximum correlation 

C,^^^(Ax)=Ci2(rMAX,Ax) 

at a particular value of the delay time tmax{Ax). Taylor's coherence hypothesis requires 
that C;(^^"^(Ax) be close to 1, and tmax{,Ax) be equal to the transit time Ax/Uq. Figure |Tl| 
shows tmax as a function of Ax/Uq, in agreement with the hne drawn for tmax = Ax/Uq 



4^ . As predicted by Taylor's hypothesis, the slope of this hne is unity. The small deviations 
are due to errors in our measurement of Ax. 

In Figure 12 we plot the maximum correlation C^'^^(Ax). As one would expect, the 



correlation decreases as Ax increases, though we have not found any simple functional 
form to flt to this decrease, nor is there to our knowledge any predicted form. The loss of 
correlation is due to the dynamic evolution of the velocity fluctuations and sets a limit to 
Taylor's hypothesis, which we quantify by deflning an "evolution length" 6^ as the separation 
for which the correlation drops to 50%. For our experiment we flnd Se — 7 cm, corresponding 
to an evolution time Xe = 6e/Uo ~ 40 ms. This length is much larger than the relevant lengths 
of the turbulent velocity fleld, as we show next. 



B. Testing Taylor's Statistical Hypothesis 

The statistical study of turbulence is framed in terms of velocity correlation functions, 
structure functions, and energy spectra; here we focus on the structure function. The 
longitudinal velocity difference between two points separated by a distance r is written as 

Sv{r, t) = (v(xi +r,t) - v(xi, t)) ■ f , 

where the unit vector f is in the downward direction of the flow. The nth order structure 
function is deflned as: 

^„(r) = ((5t;(r,t)r). 

In Figure |13| we show the structure functions 5'2(r), S^^r), and Sq^v), calculated from single 
point velocity measurements using Taylor's hypothesis: r = Uqt (solid circles). We also 
plot direct spatial measurements of these structure functions made using two probes (open 
squares). This is a direct conflrmation of Taylor's statistical hypothesis, which is one of the 
central results of our study. Note that 5*2 (r) shows a scaling region of about a decade where 
S2{r) oc r^'^, in good agreement with other experiments on turbulent soap films [^,0]. We 



also find approximately that S/^lr) oc r^'^ and SqIt) oc r^'°, as shown in the figure. The third 
moment Sr,(r) has been treated in detail elsewhere 1431 
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Some comment must be made on the observed values of the exponents, which are so 
different from S'„(r) oc r", the theoretical expectation for the enstrophy cascade in 2D 



turbulence [^,|I^. It is well known that in 3D turbulence the scaling law exponents of 
the nth order structure functions deviate from their expected value of n/3 as n gets large 
45| . This systematic difference is attributed to the intermittency of the fluctuations 0. 



However, the third order structure function must scale as 5*3 (r) ~ r even with intermittency, 
as can be derived directly from the Navier-Stokes equations |^. An equivalent derivation 
for 2D turbulence would not be relevant to the third moment in the enstrophy cascade 
range discussed here. In fact is observed to be approximately zero in this range (and 



is positive for large r) This observation suggests to us that the scaling exponents of 

the enstrophy range are more sensitive to intermittency than in 3D turbulence. Evidence of 
intermittency, indicated by non-Gaussian velocity fluctuations, has been reported previously 
for our experiment 



For r > 1 cm, the S^ir) saturate to constant values, which for S'2(r) is equal to 2t>^^g. 
This occurs roughly at the integral or outer scale 0, which characterizes the largest scales 
on which the velocity is correlated. The integral scale is defined as 

'Kr)drlvl,,s, (2) 

where h{r) is the velocity correlation function: 

h{r) = {v'{x, t)v'{x + r, t)) = 4a/s - ^'^2(r). 

At Xi = 8 cm we find io = 0.6 cm, which is much less than the evolution length 5e = 7 cm. 
Thus for the turbulence in our soap film, Taylor's hypothesis is justified: the two signals are 
correlated better than 90% for scales r < Eq (see Figs. |T^ and |T3|). 



C. Detailed Study of the Velocity Decorrelation 

There are in general two reasons for the breakdown of Taylor's hypothesis: the entrance 
of new structures into the line of travel, introducing new fluctuations into the signal, and 



the evolution of the velocity field itself. In Figure |TJ we show an overlay of the velocity 
measured at Xi = 8.0 cm vs. t, and the velocity measured at X2 = 10.0 cm vs. t — tmax 
(Ax = 2 cm). For a perfect correlation {C^^'^ = 1) the two curves would fall on top of each 
other. The arrows indicate fluctuations which have either appeared or disappeared during 
the transit time between the two probes. In effect this means that information is being 
generated, and this "new information" is partially responsible for the velocity decorrelation 
(Fig. 5). 

To explore the details of this process, we measure the coherence spectrum Cs{f) of 



the fluctuations at xi and X2 |4^. If Taylor's hypothesis of coherence were justified, then 
vi = v{xi,t) would be identical to V2 = v{x2,t — tmax)- If ^i(/) is the complex Fourier 
transform of Vi, then the standard power spectrum is Psi{f) = {vi{f)vl{f)), where v* is the 
complex conjugate of v. The coherence spectrum is 
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normalized so that Cs = 1.0 for frequencies where the two time series are coherent. A 
measurement of Cs{f) with increasing probe separation shows which modes in the turbulent 
spectrum persist longer and which evolve faster. 

Consider first the velocity power spectrum at a single point, shown as the thin line in 



Figure |I5|. This spectrum, according to the standard picture of 2D decaying turbulence 
1^,|T6|, should have a power law dependence Ps{f) ~ in the enstrophy cascade range, 
with a = 3. In an earlier soap film experiment this exponent was found to be measurably 
larger than 3; here in we find a = 3.6 ± 0.2. We compare this to the coherence spectrum, 
which we expect to be nearly equal to 1.0 for small separations. The coherence spectrum 
for Ax = 0.2 cm is also shown in Figure ^ (thick line). We see that Cs{f) is indeed close 
to unity over most of the frequency range in which the power spectrum appears. However, 
the coherence drops below Cs ~ 0.9 at / ~ 300 Hz, around the middle of the range where 
Ps{f) ~ /~", and for / ~ 850 Hz, where the power spectrum is reaching the noise floor in 
our measurement, Cs ~ 0.2. Thus already at Ax = 0.2 cm it appears that the high frequency 
components are the most rapidly evolving. In contrast, for the 2D enstrophy cascade it is 



expected that the 'eddy turnover time' is independent of size the observed falloff in 
Cs{f) may indicate viscous dissipation effects. 

The coherence spectra at five increasing separations Ax are shown in Figure |TB|, as a 
linear- log plot. In each case we see decorrelation at higher frequencies (smaller scales), while 
the low frequency part remains at a constant value which decreases as Ax increases. This 
constant correlation is approximately equal to Cf^^^{Ax), which means that the overall 
coherence of the velocity field is determined mainly by the low frequency components. The 
high frequency decorrelation also moves to lower frequencies as Ax increases. To see whether 
the whole shape of the Cs{f) follows the decay of C^"^"^(Ax), we normalize the coherence 
spectra as Cs{f)/Cf2^^ in Figure |T^. The curves do not lie entirely on top of each other, 
indicating that the cutoff at high frequencies follows a different evolution than Cil^'^ . The 
cutoff shape is well described as Cs{f) ~ log{l/ f), shown as the straight lines drawn through 
the data. This advancing cutoff is reminescent of the loss of predictability in the spectra of 



atmospheric turbulence simulations |[l^j47|| 



At larger separations (Ax > 12 cm), we find that the spectral position of this cutoff no 
longer moves to lower frequencies as Ax increases. In Figure |18| we superpose the coherence 
spectra for Ax from 12 to 22 cm, normalized by Cf^^^ (Ax). The curves lie reasonably on top 
of each other, which means that the entire coherence spectrum follows the overall decrease 
of C^"^'^ . By taking the intersection of the logarithmic fit with the line Cs{f) = Ci2^^ (Ax) 
in Figs. |1^ and |1^, we use the frequency fd of the intersection to characterize the position of 
the spectral cutoff [Q. We plot this frequency in Figure |TP|. Up to Ax ~ 12 cm, the advance 
of fd to lower frequencies is consistent with the scaling fd ~ Ax^^^^, which is slower than 
that seen for wavenumber cutoff in 2D numerical simulations of atmospheric predictability 



I^ , |47[] . For Ax > 12 cm, fd reaches a constant value of about 30 Hz, corresponding to a 



length of 6 cm, which is the size of our system (the channel width W = 6 cm). 



D. Turbulent Predictability and Taylor's Hypothesis 

The failure of Taylor's hypothesis in our experiment is closely related to the question of 
predictability in 2D turbulence pTj-^]. The general study of predictability in turbulence 
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(see for example |1T^ , |50| , [51|] ) is of particular importance to the weather prediction problem 
. Here we briefly show how our analysis parallels this general framework. Note that here 

we are comparing a developing turbulent velocity with its initial state, whereas studies of 

predictability treat the diverging evolution of two nearly identical initial states. Nonetheless 

there are several similarities between the two. 

Following Metais & Lesieur, we first define the time series of the velocity difference, or 

error time series which for our experiment is written 

Af (Ax, t) = vi{xi,t) - V2{xi + Ax, t - tmax)- (4) 

For a perfectly correlated signal this time series would be identically zero. We define the 
difference energy AE{Ax) = ((Af (Aa;, t))^), which is the second moment of Av and thus 
a kinetic energy associated with the difference series; it is analogous to the error energy in 
predictability studies [^,1^. In Figure [20| we plot a dimensionless AE, namely 



,^ , AE{Ax) , . 

^IRMS "T ^2RMS 

as a function of the decay time Ax/Uq. The function p is defined to increase from to 
1 as Ax increases, and acts as a sort of distance function between the two velocities. By 
comparing Eqs. |I| and ^ one sees that p and C12 are simply related. 

The inset to Fig. ^ shows an enlargement of p{Ax/Uo) near Ax/Uq = 0. There is no 
clear linear portion in the plot which would correspond to an exponential error increase. In 
the study of turbulent predictability an exponential error growth is used to define a sort 
of Lyapunov exponent p3[, with the error energy serving as a metric for evaluating the 



distance between co-evolving turbulent states. The data shown in Fig. ^ are in fact better 
described by a power law p ~ {Ax/UoY^'^, shown in Fig. This apparent square root 
dependence should be interpreted only as a power law dependence: the actual value of the 
exponent depends on the choice of metric function, Eq. ^ Since an exponential growth of 



the error energy depends on the linearization of an underlying equation for p, Fig. ^ may 
indicate the presence of higher order terms, similar to the quadratic saturation term used 
by Lorenz to fit error growth in an iterated map . 



The predictability time Tp is a standard measure of the time beyond which one can no 
longer project the state of a turbulent system H^JSlI]. The exact definition of such a time 



is somewhat arbitrary, though it is usually much larger than the large scale eddy turnover 
time for the turbulent flow |^. We can characterize the long time growth of p(r) by 
quantifying how long the velocity V2 remains similar to vi. The predictability time used 
by Metais & Lesieur was defined by p{Tp) = 0.5 we define an evolution time Tg such 



that p(Te) = 0.5, and find that ~ 25 ms. This is of the same order as our decorrelation 
time Te ~ 40 ms given by C^l^^ (Fig- IH)) a result which is not surprising given that the 



two functions are related. The analogy between the loss of predictability and the failure of 
Taylor's hypothesis is in fact rooted in a common cause: the loss of velocity coherence due 
to turbulence. Whether any implications can be drawn from this connection remains to be 
seen. 



V. DISCUSSION 
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A. Comparison with 3D Measurements 



We have measured the breakdown of Taylor's hypothesis for decaying turbulence in a 
flowing soap film and shown that the hypothesis is a valid assumption for statistical mea- 
surements of the turbulence (the structure functions). How do our measurements compare 
to similar experimental studies of 3D decaying turbulence? Of the six studies which to 
our knowledge provide information comparable to Fig. 12 [18-2^ we will examine three in 
detail [pT|-|23[]. Two of these studies used hot wire anemometry pl| , |22[] , and thus additional 
techniques were required to compensate for the wake of the upstream probe. Champagne 
et al. used a 'grid' made of 12 parallel channels (spacing M' = 2.54 cm) in a wind 
tunnel with a mean speed of 12 m/s, and Rcm' = 21,000. Cross-correlation measurements 



started at xi = 259 cm, where It — 0.018 and io ~ 4.2 cm. Comte-Bellot & Corrsin \T2 
made measurements behind a standard grid (M = 5.08 cm) in a wind tunnel with Rcm = 
34,000. The two hot-wire cross correlation measurements were made starting at Xi = 210 
cm, where It — 0.022 and £o — 1-1 cm. Cenedese et al. used a nonintrusive laser Doppler 
anemometer similar to our LDV (see [^), but did not use a standard grid; the turbulence 
was produced by a combination of a honeycomb and the channel walls. Their measurements 
were made in a water channel (height h = 2 cm) starting at xi = 14 cm, where It — 0.13 
and io ~ 1.0 cm. Their Reynolds numbers were also significantly lower: Rch = 4,800. 

Are there any differences between Taylor's hypothesis in our approximately two- 
dimensional soap film and in these 3D experiments? We address this question by plotting 
C^^'^ {Ax) from these three studies along with our measurements in Figure ^ The inde- 
pendent variable in this plot is Ax in units of the integral scale io- One might expect the 
decorrelation to occur more slowly in the soap film due to the absense of vortex stretching. 
However, as the turbulent intensity in our experiment is high {It = 0.14) compared to the 
two wind tunnel experiments {It ~ 0.02), our data should be directly compared only to that 
of Cenedese et al. {It = 0.13). In this case we see that indeed the correlation in our soap 
film extends to much larger values of Ax/io than in their 3D experiment. Note that Cf^^^ 
from the wind tunnel experiments also extends to much larger values of Ax/io than the 
data of Cenedese et al, probably due to the fact that their turbulent intensities are much 
lower. More definitive conclusions would come from a single experiment (in 2D or 3D) which 
measures Cf|"^"^(Ax) for several different It- 



B. Detailed Shape of the Cross Correlation Ci2(t) 

As of yet there is no rigorous underpinning to Taylor's hypothesis which would allow 
for the calculation of higher order corrections to velocity correlations, though an intriguing 



suggestion was implemented in |56]. To provide detailed information for some future theory, 
we focus on the shape of the cross-correlation function Ci2(t, Ax) around tmax- This shape 
is by definition (Eq. 0) the average convolution of a velocity fluctuation taken with itself 



a time tmax later. In Figure |2^ we show as an example Ci2(r) for Ax = 4 cm, along 
with a Gaussian distribution centered on tmax- We find that the shape is always nearly 
Gaussian, with a slight skewness around tmax consistently towards the positive. The width 
of the Gaussian does not broaden as Ax increases, though the maximum does decrease as 



shown in Figure |2J. Thus the development of the cross-correlation cannot be treated as a 
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diffusion-like process, for which the width would increase as the maximum decreases. The 
small positive skewness is also not strongly dependent on Ax. 

VI. CONCLUSION 

In this paper we focused on the breakdown of Taylor's coherence hypothesis in a turbulent 
soap film, a quasi-2D experimental system. We have shown that for the lower order moments 
the statistical hypothesis works well, even when the actual cross correlation between the 
two probes is low. As the relevant length scale of this decorrelation is much larger than the 
integral scale of the turbulence (5e >> f-o), this phenomenon is outside the region usually 
considered by most studies: it is the turbulence beyond the scaling range. Yet this evolution 
contains untapped information, as we have indicated. The failure of Taylor's hypothesis 
may thus shed light on deeper problems in turbulence. 
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FIGURES 



FIG. 1. A diagram of the experimental setup: A) upper soap reservoir; B) hooks which hold 
plastic wires between which soap film flows; C) comb behind which turbulence is generated in the 
film; D) lower soap reservoir. The two LDV probes which measure the velocity in the film are 
labeled 1 and 2. Downstream distances are labeled as referred to in the text. 



FIG. 2. The decay of the turbulent intensity It behind the comb in our soap film: a) It vs. down- 
stream distance x; b) the same data plotted as It~^'^ vs. x/M. The straight line is a linear fit (see 
text). 



FIG. 3. Simultaneous traces of the velocity vs. time at xi = 8.0 cm and X2 = 8.5 cm behind 
the comb. 



FIG. 4. The cross correlation Ci2(t, Ax) vs. delay time r for several different values of Ax (as 
labeled) . 



FIG. 5. The delay time tmax for the maximum of the cross correlations in Fig. |10| vs. the 
transit time Ax/Uq. The solid line corresponds to tmax = Ax/Uq. 



FIG. 6. The maximum value of the cross correlations in Fig. |T^, 0^2^-^ vs. the downstream 
separation Ax. The arrow shows the evolution length 5e~^ cm (see text). 



FIG. 7. A diagram of the experimental setup: A) upper soap reservoir; B) hooks which hold 
plastic wires between which soap film flows; C) comb behind which turbulence is generated in the 
film; D) lower soap reservoir. The two LDV probes which measure the velocity in the film are 
labeled 1 and 2. Downstream distances are labeled as referred to in the text. 



FIG. 8. The decay of the turbulent intensity It behind the comb in our soap film: a) It vs. down- 
stream distance x; b) the same data plotted as It^'"^ vs. x/M. The straight line is a linear fit (see 
text). 
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FIG. 9. Simultaneous traces of the velocity vs. time at xi = 8.0 cm and X2 = 8.5 cm behind 
the comb. 



FIG. 10. The cross correlation Ci2{t, Ax) vs. delay time r for several different values of Ax 
(as labeled). 



FIG. 11. The delay time tmax for the maximum of the cross correlations in Fig. |l^ vs. the 
transit time Ax/Uq. The solid line corresponds to tmax = Ax/Uq. 



FIG. 12. The maximum value of the cross correlations in Fig. 10, 0^2^^ vs. the downstream 
separation Ax. The arrow shows the evolution length 5e — ^ cm (see text). 



FIG. 13. Experimental check of Taylor's "statistical hypothesis": a) second order structure 
function 5*2 (r) taken at x = 8 cm using Taylor's hypothesis; b) fourth order structure function 
Si{r); c) sixth order structure function 5*6 (r). The open squares are direct spatial measurements 
made using two probes, and the lines correspond to the fitted scalings described in the text. 



FIG. 14. Overlay of the velocity measured at xi = 8.0 cm vs. t, and the velocity measured at 
xi = 10.0 cm vs. t — Tmax- The arrows indicate fluctuation spikes which have either appeared or 
disappeared during the transit between the two probes. 



FIG. 15. A comparison of the coherence spectrum (Eq. ^ for Ax = 0.2 cm, and the power 
spectrum taken at xi = 8 cm. The line is a fit to the power law Psi{f ) ^ f~" with a = 3.6 ± 0.2. 



FIG. 16. Coherence spectra for several values of Ax as a linear-log plot of /. 



FIG. 17. Coherence spectra as in Fig. It normalized by the maximum cross correlation C^"^^ 



from Fig. 12. The straight lines through the high frequency part of the coherence spectra corre- 



sponds to a logarithmic decay law Cs{f) ~ log{\/ f). 
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FIG. 18. Coherence spectra for several large values of Ax as a linear-log plot of /, normalized 
by the maximum cross correlation Ci2^^ from Fig. 12. 



FIG. 19. The frequency scale fd characterizing the logarithmic decay Cs{f) ~ log{l/f) (lines 
in Fig. 17) vs. the separation Ax. The straight line represents the scaling law ~ Ax~^/^. 



FIG. 20. The normalized error energy p{Ax/Uq) as a function of time Ax/Uq past the reference 
point xi, shown as a log-linear plot. The arrow indicates the evolution time Te as defined in the 
text. Inset: an expanded view near the origin. 



FIG. 21. A log-log plot of p{Ax/Uq) as a function of downstream time Ax/Uq. The straight 
line corresponds to a power law dependence (Ax/C/q)^^^ (see the text). 



FIG. 22. A replot of C^|^^ from Fig. p as a function of probe separation Ax normalized 



by the integral scale Iq. 
Ref. |3| 



Also shown for comparison are 3D results from Ref. |^], Ref. |22], and 



FIG. 23. Gaussian fit to the shape of Ci2(r) at Ax = 4 cm. 



FIG. 24. An overlay of several Ci2(t) at various Ax (as shown), plotted vs. r — tm AX- 
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